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Abstract 

We study SUSY phenomenology of the KKLT (Kachru-Kallosh-Linde-Trivedi) 
type scenarios of string theory compactifications with fluxes. This setup leads to a 
specific pattern of soft masses and distinct phenomenological properties. In partic- 
ular, it avoids the cosmological gravitino/moduli problems. Remarkably, the model 
allows for the correct abundance of SUSY dark matter consistently with all experi- 
mental constraints including the bound on the Higgs mass, b — > sj, etc. This occurs 
for both small and large tan /?, and requires the SUSY spectrum above 1 TeV. 
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1 Introduction 



String compactifications with fluxes have recently attracted considerable attention. 
The presence of fluxes allows to stabilize most moduli and eliminate these unwanted 
scalars from the low energy action [1]. One of the most attractive setups in which 
all the moduli are fixed and the cosmological constant is zero or small is a model 
due to Kachru, Kallosh, Lindc and Trivcdi (KKLT) [2]. The consequent SUSY 
spectrum exhibits a number of interesting features [3,4]. In particular, the soft terms 
receive comparable contributions from gravity (modulus) mediated [5] and anomaly 
mediated [6] SUSY breaking^. Another robust feature is a hierarchy among the 
MSSM soft masses, the gravitino and moduli masses, 

"^MSSM < "^3/2 < "Imoduli • (1) 

Some phenomenological aspects of this class of models have recently been studied 
in Refs. [9,10]. In particular, it was observed that the heavy gravitino and moduli 
alleviate cosmological problems associated with late decays of these particles [10]. 
Also, the pattern of soft masses was found to be quite distinct [9]. 

In the present work, we undertake a comprehensive study of phenomenological 
properties of the model. We analyze experimental constraints on the spectrum from 
collider bounds on sparticle and Higgs masses, BR(6 57), etc. as well those im- 
posed by correct electroweak symmetry breaking and absence of charge and color 
breaking minima in the scalar potential. Then we study compatibility of these con- 
straints with the requirement of the correct SUSY dark matter abundance. Although 
the spectrum is very constrained and parametrized in terms of three continuous quan- 
tities only (m3/2, a and tan/3), we find that the right amount of dark matter can 
be produced in considerable regions of parameter space. Unlike in the common 
mSUGRA model, both low and high values of tan/3 are allowed. 

The outline of the paper is as follows. In section 2 we introduce the KKLT model, 
in section 3 we analyze the consequent soft SUSY breaking terms and in section 4 
we study relevant phenomenological constraints. Our conclusions are presented in 
section 5. Some technical details concerning the anomaly mediated soft terms are 
given in the Appendix. 

2 The KKLT setup 

In this section we discuss the KKLT construction and its main features. The KKLT 
setup is based on Calabi-Yau compactifications of type IIB string theory with fiuxes 
[11]. The presence of background fluxes in the compactified space, that is non-zero 
vacuum expectation values of certain field strengths, allows one to fix all complex 
structure moduli as well as the dilaton [1]. The former parametrize the shape of the 
internal manifold and in the absence of fiuxes have a zero potential to all orders in 
perturbation theory. Internal fiuxes create a potential for moduli thereby mitigating 
a number of phenomenological problems associated with light or massless moduli. 

similar pattern also appears in the heterotic string [7]. For phenomenology of compactifications 
with fluxes, see also [8]. 
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This mechanism, however, does not apply to the overall T-modulus parametrizing 
the size of the compact manifold. The KKLT proposal is to invoke nonperturbative 
mechanisms such as gaugino condensation on D7 branes to stabilize the remaining 
modulus. As a result, the vacuum energy in such a theory is negative which requires 
further modifications of the setup. To this end, KKLT add a contribution from a 
non-supersymmetric object (anti-brane) which does not significantly affect moduli 
stabilization. Thus the setup requires the presence of a number of D7/D3 branes 
and an anti D3 brane. The final outcome is that (i) all moduli are fixed, (ii) the 
cosmological constant is small and positive. This is the major achievement of the 
model. 

Let us now consider the KKLT model in more detail. We start with a 4D super- 
gravity scalar potential. A supergravity model is defined in terms of three functions: 
the Kahler potential the superpotential W ^ and the gauge kinetic function /. The 
scalar potential is given by 

^SUGRA = [k'^DiWDjW* - 3|Vr|2) . (2) 

Here DjW = djW + WdjK is the Kahler covariant derivative of the superpotential 
and K^'^ = (djdjK)^^ . The gravitino mass is given by 

m3/2 = (3) 

and the SUSY breaking F-terms are 

= -Mpi^ e^''^K^^DjW* . (4) 

Given K and W as functions of the fields in the system, one minimizes the scalar po- 
tential VsuGRA and finds whether supersymmetry is broken {Fj 0) in the vacuum 
or not. In supergravity, vanishing of the cosmological constant imposes the relation 
rriy^ ~ K^'^ FjFj, therefore 171^/2 serves as a measure of SUSY breaking. The grav- 
itino acquires its mass through the super-Higgs effect, that is, it absorbs the spin 
1/2 Goldstino associated with spontaneous SUSY breaking. The MSSM soft masses 
are controlled by the F-terms such that typically one expects the soft masses to be 
of the order of the gravitino mass. The moduli masses are found from derivatives of 
VsuGRA at the minimum and are also within one- two orders of magnitude from 777.3/2 
(cf. [12]). 

In the KKLT setup, the total scalar potential is given by the sum 

V = FsUGRA + Mift , (5) 

where Viift is an explicitly SUSY breaking contribution which serves to lift the min- 
imum of the potential to a Minkowski or de Sitter vacuum. With a general Vjift, 
the gravitino mass Q is not related to the F-terms. and is an explicit mass term. 
Similarly, the moduli masses found by differentiating V are not related to or 
the F-terms. This has its advantages since the gravitino and the moduli can be made 
heavy so as to avoid cosmological problems associated with late decays of these par- 
ticles. At the same time, the F-terms can be kept small enough to produce a TeV 
MSSM spectrum required by naturalness in the Higgs sector. 
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Let us now consider the specifics of the KKLT scenario. VsuGRA is a function 
of the T-modulus (as well as the MSSM fields which we suppress) with the Kahler 
potential and the superpotential given by 

K = ~3ln{T + f ) , = u-o - Ce""^ . (6) 

Here T is related to the compactification radius R, RcT ~ R^, wq is a constant 
induced by the fluxes, C is a model-dependent coefficient and a is related to the beta 
function of gaugino condensation on the D7 branes, a = Svr^/iVc for SU{Nc). The 
lifting potential due to the presence of the anti D3 brane is 

^lift = J^q:fy^ ' (7) 

with n being an integer (n = 2 in the original KKLT version) and D is a tuning 
constant allowing to obtain a Minkowski/de Sitter vacuum. 

At D = 0, the scalar potential is minimized at V = —Stji^^^^pi- The addition 
of the lifting term leaves the value of T at the minimum essentially intact. This 
is because the supergravity potential is exponentially steep unlike the lifting term. 
Thus, the effect of the lifting term is simply to change the vacuum energy to a 
small positive or zero value. This is achieved with D ~ ml^^M^^ ~ lO^^'^M^j. 
Such a small value may appear unnatural. However, one should remember that 
the background geometry in the KKLT model is warped, ds^ = e^"^^y\]fj,iydx^dx'^ + ... 
with y parametrizing the compact dimensions and A{y) being a flux-dependent warp 
factor. At the location of the SM fields the warping can be negligible, e^^ ~ 1, 
whereas at the location of the anti D3 brane the space can be significantly warped, 
^2A <^ i_ In ^]2is case, the natural mass scale on the anti D3 brane is much smaller 
than the Planck scale and can be chosen to be of the order of the intermediate scale, 
m ~ e'^Mpi ~ y^rn^j2Mp\. Thus the desired value of D can be obtained by placing 
the anti-brane at the appropriate point in the compact space. 

Minimizing the scalar potential one finds, 

-3/2 - Mp, ^^^^ , 

a ReT ~ - In (7713/2 /Mpi) , 

Ft "^3/2 
ReT ^ a ReT ' 
niT ~ a ReT 771,3/2 • (8) 

To get a TeV MSSM spectrum, wq should be adjusted to be very small, lO^"*^^, 
which can be achieved by finetuning fiuxes in the underlying string theory. Then, 
a ReT ~ 25. This is a moderately large parameter leading to a hierarchy among the 
gravitino, the modulus and the MSSM soft masses. Indeed, assuming that the MSSM 
fields live on D7 branes, the soft masses arc contr oiled by F'^Ot Ini^MSSM ~ F'^ /ReT 
[13] and thus are suppressed by a ReT compared to the gravitino mass. On the other 
hand, the modulus mass is enhanced by the same factor compared to 7713/2- As 
mentioned earlier, this moderate hierarchy is highly desired from the cosmological 
perspective: the modulus and the gravitino produced in the early Universe would 
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decay before the nucleosynthesis and thus would not affect the abundances of light 
elements. 

We note here that the most important effect of the presence of an anti-brane is 
lifting the vacuum energy. As argued in Ref. [4], other effects due to the existence of 
an explicit gravitino mass term on the anti D3 brane or other explicit SUSY breaking 
terms are expected to be suppressed by warping. 

Concerning localization of the MSSM fields, there are a few options: they can live 
on D7 branes, D3 brancs or on D7 and D3 branes^. There are certain advantages 
and disadvantages to each of these choices. If the observable fields are localized on 
the D3 branes, the MSSM spectrum is plagued by negative slepton masses squared - 
the usual problem of anomaly mediated SUSY breaking. Furthermore, it is difficult 
to get (semi-) realistic quark/lcpton flavour structures. In the case of D7 brancs, 
such problems do not arise: SUSY breaking is communicated by both the anomaly 
and the modulus F-term such that all masses can be made positive. For the flavour 
structures, in principle one can use the successful technology of intersecting branes 
[14]. On the other hand, the theoretical calculations are not well under control since 

ReT = ~ 2 (9) 

S'gut 

requires a non-perturbative string coupling (cf. Eqs.(2),(6) of Ref. [4]). In any case, 
there are still some outstanding theoretical issues in this setup which have to do, for 
instance, with effects of explicit SUSY breaking contributions. We will not attempt 
to resolve these problems here. Instead, we will use the KKLT scenario as motivation 
to study certain patterns of soft SUSY breaking terms. As argued above, the setup 
with the MSSM on D7 branes is phenomenologically more appealing and we will take 
it as an assumption. 



3 The MSSM soft terms 

In this section we discuss the soft supersymmetry breaking terms for the MSSM fields 
living on D7 branes in the KKLT setup. Their main feature is that they interpolate 
between the soft terms of anomaly-mediated SUSY breaking^ and those of gravity- 
mediated SUSY breaking. This pattern appears generically whenever moduli are 
stabilized close to a supersymmetric point and leads to distinct phenomenology. 

The Kahler potential and the kinetic function for the MSSM gauge fields are given 
by [13] 

i^ = -3in(r + T) + ^^J|t- , /„ = r. (10) 

Here Qi are the MSSM matter fields and a = 1,2, 3 runs over the GUT normalized 
C/(l), SU{2) and SU{2)) group factors, = {0, 1/2, 1} are constants depending on 
the origin and localization of the matter fields. For definiteness, we will set = 

^In the original KKLT proposal, the MSSM fields were implicitly assumed to be localized on D3 branes. 
As we mention below, this choice is problematic due to the presence of tachyonic sleptons. 

^The anomaly-mediated contribution is usually present in string models [15], but may be absent in 
certain cases [16]. 
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Figure 1: The KKLT setup with the MSSM fields on D7 branes. The factor e^^^y^i represents warping 
along the compact dimension y (Klebanov-Strassler throat) . 



in what follows. The choice Ui = 1 would lead to tachyons, whereas phenomenology 
of the Ui = 1/2 model would be quite similar to that of the nj = case (with larger 
tachyonic areas in parameter space). 

The gaugino and soft scalar masses are generated by the auxiliary component Ft 
of the modulus superfield. Their magnitude is controlled by Ft ~ ?7Z3/2/a <SC "^3/2- 
These tree level terms are much smaller than the gravitino mass and are comparable 
to the loop-suppressed anomaly-mediated contributions. The scale of the anomaly- 
mediated contributions is set by F$/167r^, where F<j> = 171^/2 + ^FtBtK ~ ?7^3/2 
and $ is the conformal compensator. It is convenient to parametrize our F-terms in 
terms of a new scale Mg defined by 

F$ ~ 7713/2 = Wtt^Ms , =aMs. (11) 

Here a depends on the shape of the lifting potential and is given hy a ^ IQn^ 2aReT 
(note the difference from a defined by Choi et al. in Ref. [9]!). Its precise value 
depends on further details of the model such as the string scale, the gravitino mass, 
etc^. For the original KKLT lifting potential n = 2 and a lies in the range 4.8 -r- 6. 
With other choices of the lifting potential, different values of a can be obtained, e.g. 
a = 7 -\r 9 for n = 3. In the limit a — > we recover pure anomaly mediation, while 
a S> 5 corresponds to gravity (modulus) mediation. 

The soft terms in the mixed anomaly-modulus mediation scenario are controlled 

^This dependence appears since aReT which solves the equation e^"^°^ = 2aRc t Ij" depends on C 
and the gravitino mass. 



5 



by the scale Mg and given by [4]: 



Ma = M,[a + bagl] 



m: 



2 - M^[a^ -ji + 2aiT + T)dTli 



Aijk = M, [3a - 7i - 7j - 7fc] + A^ij-fc . (12) 

Here ba are the beta function coefficients for the gauge couplings ga, 7i is the anoma- 
lous dimension and 7j = STr^ g^^^ . In supersymmetric models, 

7. = 2Y,glci{Q^)-Y,\y^\' ^ 

a Hi 

where C2{Qi) is the quadratic Casimir corresponding to the gauge representation 
of Qi. In the second term, the sum runs over all physical Yukawa couplings yi 
involving Qj. The coefficient by. describes the running of the Yukawa couplings, 
= j^Viby^. Finally, the scalar soft masses contain a mixed anomaly-modulus 
contribution proportional to dxji which appears due to the T-dependence of the 
gauge and physical Yukawa couplings, 

{T + T)dTli = -2^glC^{Qi) + 3^ [y^p + Si . (14) 

a i 

All relevant RG parameters are listed in Appendix A. Numerically, the anomaly and 
gravity mediated pieces in Eq. H12|) are roughly the same at a < 3. 

6i in Eq. (|14|) and AAjj^ in the expression for the A-terms account for a po- 
tential T-dependence of the Yukawa couplings^, A^y^ oc STlnl^j^. The presence 
of this term as well as its specific form depend on the theory of flavour and can- 
not be analyzed in full generality. For simplicity, we will omit this term in most of 
our phenomenological analyses, yet we will comment on some of the effects it can 
generate. 

The remaining two parameters important for SUSY phenomenology are the /x and 
the Bn terms. Since these are responsible for electroweak symmetry breaking, their 
magnitude is bounded by the scale of the soft masses, that is around 1 TeV. This is 
rather difficult to achieve in models similar to the anomaly mediation scenario since 
the natural value for the i?-term would be ~ "^3/2 ^ ^s- Nevertheless, the 
desired values can be obtained with some finetuning given the //-term is generated in 
the superpotential, AW = KH1H2, as well as the Kahler potential, AK = k' H1H2 + 
h.c. Then, parametrizing /i as 



= l^w + IJ'K , (15) 

one has 

= cim3/2Hw + C2ms/2fJ-K , (16) 

with ci_2 being order one constants which depend on k and k' . Adjusting and 
appropriately, /i and B of order Ms can be obtained [9] . The practical conclusion 



'In simple cases, the holomorphic Yukawa couplings are T-independent [17]. 
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is that, lacking a compelling model of generating /x and B, they should be treated as 
adjustable parameters so as to produce correct electroweak symmetry breaking. 
Let us now overview main features of the resulting SUSY spectrum. 

(i) Moduli/gravitino problem. A characteristic feature of the spectrum is a 
moderate hierarchy (a factor of 30 or so) between the MSSM soft masses and the 
gravitino mass as well as between the gravitino mass and the moduli masses. As 
discussed in Ref. [10], this is advantageous from the cosmological perspective since 
the gravitino and moduli are heavy enough to decay before the nucleosynthesis and 
not to affect abundances of light elements. 

(ii) Tachyons. Pure anomaly mediation is notorious for its negative slepton 
mass squared problem. In the KKLT setup, there is an additional gravity mediated 
contribution which rectifies the problem. The absence of tachyons imposes a lower 
bound on the parameter a. Indeed, the GUT scale boundary condition for the slepton 
masses of the first two generations reads 

ml « [-1 - 2a + a^) , 

ml « (-2 - a + a^) . (17) 

To avoid tachyonic sleptons, a > 2 is required. For the squarks, 

mQ « (2 - 4a + a^) M"^ , 
mfj « (l - 3a + a^) M| , 

ml « (2 - 3a + a^) . (18) 

Although the squark masses are positive in pure anomaly mediation, they become 
tachyonic^ for 0.5 < a < 4 due to the mixed anomaly-modulus contribution propor- 
tional to a. In conclusion, the tachyons which signify color or charge breaking minima 
are absent for a > 4. This bound has important implications for phenomenology. In 
particular, most of the parameter space with characteristic signals of anomaly me- 
diation such as a wino LSP is excluded. Curiously, a ~ 5 predicted by the original 
KKLT model is on the safe side. 

(iii) LSP. In the non-tachyonic region, the bino is the lightest gaugino. Our 
numerical analysis shows however that for 4 ^ a ^ 8 the LSP is dominated by the 
Higgsino component. This can be explained as follows. The anomaly and modulus 
mediated contributions add up in the GUT scale bino mass. Mi Ms{a + 33/10) 
but partially cancel in the gluino mass M3 Ri Ms{a — 3/2). It is well known that 
the low energy value of the Higgs mass parameter ?n.|^^ and, consequently, the fi- 
term, is typically controlled by the GUT scale gluino mass, //^ « (2 3)M^. Thus, 
for intermediate a where the suppression of the gluino mass is effective, we get 
1^1 < Mi(TeV) ~ 0.4Mi and the LSP is higgsino-like. This is certainly desired from 
the SUSY dark matter perspective. We also note that the stau LSP is not possible 
in this scenario since unlike in mSUGRA the scalar and the gaugino masses cannot 
be varied independently. 

®For a > 2, the squark masses squared are positive at the EW scale due to the RG running. However, 
2 < a < 4 lead to tachyonic squarks at the GUT scale which signifies existence of color breaking minima 
in the effective potential. 
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(iv) Mirage unification. An interesting feature of the scenario is the occurence 
of mirage unification [9]. That is, even though the gaugino and the scalar masses 
do not unify at the GUT scale, RG running of these quantities makes them unify at 
some intermediate scale. Indeed, the solutions to the 1-loop RG equations (neglecting 
Yukawa contributions) read 



Mail^) = Ms 



1 _ ^<^9qut 1q„ _ 



GUT 



m?(^) = M,V(^l + 2^2M)-2^^Mf(^). (19) 



M^.Te-'^'/" , (20) 



At the mirage scale fJ-mh, 

/^mir — ^^^GUT 

all gaugino and scalar masses of the first two generations unify, 

M^{^^^ir) = mji^i^ir) = {aMsf . (21) 

This is truly a mirage scale as there is no physical threshold associated with it. 
Furthermore, the third generation scalar and the Higgs mass parameters do not 
unify at that scale. We note that for a ~ 5 the mirage unification occurs at an 
intermediate scale, /imir ~ 10^^ GeV. In this case, the low energy spectroscopy is in 
some respects similar to that of gravity mediation with an intermediate string scale. 
In particular, the hierarchy between the squark and slepton masses is reduced, as 
compared to mSUGRA models. 

(v) FCNC problem. The FCNC problem can only be addressed after realistic 
Yukawa flavour structures have been obtained. The problem appears when the rii 
parameters of Ea. HlU|) are generation-dependent, rii are generally correlated with the 
Yukawa structures [18] such that the problem might actually be absent in realistic 
models. In our analysis, we simply assume that all nj = in which case the FCNC 
are suppressed. 

In any case, as we argue in the next section, consistency with accelerator con- 
straints requires a heavy SUSY spectrum, 1-5 TeV. Since the scalar mass matrix is 
diagonal, even generation-dependent choices for rii would not lead to any significant 
FCNC problem [19]. 

(vi) CP problem. The equations of motion require Arg(F$) ~Arg(i<V). As a 
result, the CP phase in the gaugino masses is aligned with the universal CP phase 
of the A-terms. This means that the physical phases Aig{M*A) vanish. Yet, there 
remain two sources of dangerous physical phases. First, it is a phase of the type 
Arg(M*i3)= Arg^M*(S/i)/x*^. From Eas. (|15ll6|) it is clear that this phase is pro- 
portional to 

It can be associated with the phase of the //-term since Ma and B/j, can be made 
real by ?7(1)_r and C/(1)pq rotations. This phase is strongly constrained by EDM 
experiments, cpfj, < (see e.g. [20]). Since there is no reason for /i^y and fix to 

be aligned and the presence of both is required by correct electroweak symmetry 
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breaking, EDMs are overproduced unless the SUSY spectrum is heavy. We note that 
the same problem appears in the well known dilaton-domination scenario. 

The second source of EDMs is associated with A-term non-universality [21], 
namely the term AAj^jt in Ea. (|12jl . Even if the A-terms could be made real by 
a U{1)r rotation, they would be flavour-dependent. That is, they would not be 
aligned with the Yukawa matrices. The latter are necessarily complex and require 
diagonalization involving complex rotation matrices. Specifically, defining Aij = 
AijYij with i,j being the flavour indices, the Yukawas and the A-terms transform 
under a basis change as 

y^vIyvr, 

A^vIAVr, (22) 

where Vl,/? diagonalize the Yukawa matrices in the up- and down-sectors. In the 
basis where the Yukawa matrices are diagonal, the A-terms have a general form and 
their diagonal entries involve CP phases. The resulting EDMs usually exceed the 
experimental bounds by orders of magnitude [21]. In our analysis, we will simply 
assume that the dangerous term AAiji^ is absent. 

We find however that although the CP phases are present generically, the induced 
EDMs are suppressed due to the heavy SUSY spectrum (1-5 TeV) such that no 
significant CP problem exists. In what follows, we set the CP phases to zero for 
simplicity. 

To conclude this section, we find that the KKLT setup leads to an interesting 
pattern of the soft masses. Although it may not solve all the problems, it has a 
number of positive features, in particular with regard to cosmology. In the next 
section, we present our detailed numerical study of the spectrum and low energy 
observables. 

4 Phenomenology 

As discussed in the previous section, the model contains four free parameters at the 
GUT scale: the gravitino mass 171^/2 = ISvr^M^, the modulus to anomaly mediation 
ratio parametrized by a, the /i-term and the B term. The absolute value of /i is 
determined by requiring correct electroweak symmetry breaking, whereas its sign 
remains free. Further, it is conventional to trade B for a low energy parameter 
tan/3 = {H2)/{Hf), which is a function of B and other GUT scale parameters. 
Thus, the parameter space for phenomenological studies is 

m3/2 , a , tan/3 , sgn(/i) . (23) 

These are our input parameters at the GUT scale, ~ 2 x 10^^ GeV. We assume that 
effective field theory is valid below this scale and use RG equations to derive the 
low energy SUSY spectrum. This is really an assumption since the string coupling 
is large in the regime considered and the effective field theory approach may not be 
valid. To this end, we use the bottom-up perspective and study the pattern of the 
soft terms hinted by the KKLT model. 
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Once tan/3 and sgn(/Lt) have been fixed, we scan over the gravitino mass < 
"^3/2 < 150 TeV and < a < 10. The low energy mass spectrum is calculated 
using the Fortran package SUSPECT [22] and its routines described in detail in Ref. 
[23]. Evaluation of the b ^ S'j branching ratio, the anomalous magnetic moment 
of the muon and the relic neutralino density is carried out using the routines of 
micrOMEGAs 1.3.1 [24,25]. 

In what follows, we divide constraints on the model into two classes which we 
call "theoretical" and "accelerator" . The theoretical constraints include correct elec- 
troweak symmetry breaking, absence of color and charge breaking minima, as well 
as the dark matter abundance consistent with the WMAP limits. The accelerator 
constraints include bounds on the Higgs and sparticle masses, the b ^ branching 
ratio and similar observables. 

4.1 Theoretical constraints 
4.1.1 ElectrowecLk symmetry breeiking 

Minimizing the MSSM Higgs potential leads to the standard relation 

„ —m% tan^/3 + m?r 1 „ 

imposed at the SUSY breaking scale defined by the average stop mass, MsusY = 
^JrniTrTi^. In most cases, it is well approximated by 

IJ? ^ -ml^ - . (25) 

When the right hand side is negative, electroweak breaking cannot occur, ^"h,^ 

at MsusY is computed by using its RG evolution from the GUT scale, -q\^^ ~ 
(iytimjf^ + mfj^ + rnq^ + Af) with /j, being the scale parameter. The result depends 
most sensitively on the gluino mass M3 at the GUT scale which increases m^^ and 
m^g. Typically, one finds mjf^{MsvsY) ^ -{2^ 3)M|. 

In the model under consideration, the anomaly and the gravity contributions 
appear in M3 with opposite signs, M3 fs Ms{a — 3/2). For low a, the effect of M3 on 
mfj^ (Msusy) is suppressed such that other RG contributions become more important 
and a negative 'm'jj^{Ms\jsY) cannot be obtained. Thus, the requirement of correct 
electroweak symmetry breaking imposes a lower bound on a. Taken together with 
the constraint from the absence of tachyons, this bounds requires typically a > 4-=-6. 



4.1.2 Colour and cheirge breeiking minima 

Generically, supersymmetric models have many flat directions in the field space. 
SUSY breaking terms usually lift these directions, but may also induce global or deep 
minima which break the electric charge and colour symmetries (CCB minima) [26]. 
It is important to verify that such minima do not develop. 

Some of the dangerous CCB minima appear along D-flat directions when the 
trilinear A-terms are sufficiently large. Absence of such minima imposes constraints 
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on the magnitude of the A-terms. In particular, 

At<^{ml,+ml+ml) . (26) 

Eq. (|12jl imphes that this constraint is usually respected. We have also checked this 
statement numerically. 

Another type of constraints comes from F— and D— flat directions. Among the 
dangerous flat directions are those corresponding to the gauge invariants LH2 and 
LLE, LQD. Absence of CCB minima along these directions usually guarantees their 
absence along the remaining directions (see e.g. [27]). A CCB minimum develops for 
a negative rn'jj^ +m\ due to the negative and large in magnitude rn^u^ at low energies. 
We find however that this does not occur in viable regions of the parameter space, 
mainly due to the negative anomaly mediated contribution to M3 which reduces the 
magnitude of fn^u^ ■ Altogether, absence of CCB minima does not constrain the model 
significantly. 



4.1.3 Neutralino dark matter 

The 2(T WMAP limit on the dark matter relic abundance is [28] 

0.094 < VLDKih^ < 0.129 . (27) 

In SUSY models, the typical dark matter candidate is the lightest neutralino and it is 
the case here. In most of the parameter space, the lightest neutralino Xi is the LSP. 
Assuming i2-parity conservation it is stable. Then to get the consistent dark matter 
abundance one has to make sure that the neutralinos annihilate efficiently enough to 
satisfy the bound ()27() . In this computation we will assume that the LSP abundance 
is thermal. Further, we will treat regions of the parameter space violating the upper 
bound in as "ruled out", those within the bounds as "favoured" and those 
below the lower bound as "allowed". The last case implies that there are additional 
ingredients to dark matter, beyond the MSSM, or that dark matter production is 
non-thermal. 

The four neutralinos x'i=i 234 superpositions of the neutral fermionic partners 
of the electroweak gauge bosons and W3 , and the superpartners of the neutral 
Higgs bosons H^, H^. In the basis {B^, W3 , H^, H^), the neutralino mass matrix is 
given by 

—m^ cos /3 sin 6*^1/ mz sin PsinO^y \ 

M2 mz cos P cos 6^^ —mz sin P cos 9^, 

mz cos P cos 9y^, —fi 

-mz sin P cos 6^, / 

This is diagonalized by an orthogonal matrix Z such that the lightest neutralino is 
given by 

X? = ^11^0 + ZuW^ + ZisH^a + ZuH^ . (28) 

Xi is usually called "gaugino-like" if P = jZnp + |Zi2p > 0.9, "Higgsino-hke" if 
P < 0.1, and "mixed" otherwise. 



M 



N 



( Ml 


—mz cos P sin 6^, 
\ mz sin P sin 9^ 
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It is instructive first to recall the situation with dark matter in the minimal 
supergravity model (mSUGRA). In most of the parameter space, the lightest neu- 
tralino is mainly the bino and, as a consequence, the annihilation cross section is 
small producing too large relic abundance. Nevertheless, there are three corridors 
in the parameter space where the cross section is enhanced. First, there is the 
stau-neutralino coannihilation branch, i.e. the region where the stau mass is al- 
most degenerate with that of the LSP. Second, there is the A-pole region where 
4(mj^o)^ ~ ~ m'jj^ — rnj^^^ — and the dominant neutralino annihilation pro- 
cess is due to the s-channel pseudo-scalar Higgs exchange. Finally, the annihilation 
cross section is enhanced if the LSP is of the Higgsino type, which occurs for small 
/X. In that case the neutralinos annihilate efficiently through the Z boson exchange 
and also coannihilate with the charginos. 

We find that the first option cannot be realized in the model under consideration. 
The reason is that the stau is always much heavier than the LSP since, unlike in 
mSUGRA, the gaugino and the scalar masses cannot be varied independently. How- 
ever, the A-pole and the Higgsino LSP corridors are indeed present and the WMAP 
bounds are respected in considerable regions of the parameter space. 

4.2 Accelerator constraints 

4.2.1 Direct search constraints 

An important constraint on the parameters of the model comes from lower bounds 
on the sparticle and Higgs masses due to direct collider searches. We implement 
these bounds by first ensuring absence of tachyons in the squark and slepton sectors 
and then applying the LEP2 constraints. The most restrictive bounds are due the 
chargino mass constraint, > 103.5 GeV, and, particularly, due to the lightest 

Higgs mass constraint. In the decoupling limit Ma 2> Mz which is applicable in 
all of the viable parameter space, the latter bound is rrih > 114 GeV at 3a. It is 
well known that this bound is sensitive to the value of the top mass. In most of 
our analysis, we have used the central value mt = 178 GeV. We have subsequently 
studied sensitivity of the results to the precise value of the top mass by considering 
the 2(7 limiting cases, mt = 174 GeV and mt = 182 GeV. 

4.2.2 BR(b S7) 

The supersymmetric spectrum is constrained indirectly by the branching ratio of 
the b decay. The most important SUSY contributions involve the chargino- 

stop loops as well as the top-charged Higgs loops. We impose the 3a bound from 
CLEO [29] and BELLE [30], 2.33 x 10"^ < BR{b sj) < 4.15 x 10"^. We find that, 
typically, the 6 — > 57 bound is more important for /x < 0, but can also be relevant 
for /X > 0, particularly at large tan/3. 

4.2.3 Muon g-2 

The 2.7a deviation of the experimental value of the muon anomalous magnetic mo- 
ment [31] from the SM prediction [32] may be interpreted as indirect evidence for 
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physics beyond the Standard Model and, in particular, supersymmetry. This devia- 
tion favours a relatively light SUSY spectrum and a specific set of SUSY parameters, 
e.g. a positive We find however that consistency with other data requires a rather 
heavy spectrum in our model such that the muon g — 2 deviation cannot be explained 
(unless tan/? is large). Thus we will treat g — 2 simply as a 3a constraint on the 
model and will display the 2a bands where relevant. 

The discrepancy 5a = a^^^ — is measured to be (5a^ = (27.1 ± 10) x 10^ ^'^ 
if e~^e~ annihilation data for the calculation of a^J^ are used. When the tau data 
are used instead, a smaller discrepancy is found. In this case, the 3a bound is 
6a^ > -6 X 10-1°, which we use in our analysis as a bound on the SUSY contribution 
(Ja^^^'^. For /i < 0, this excludes part of the parameter space with a relatively light 
spectrum. For > 0, it imposes no constraint. In that case, we display in our figures 
the 2a band 5af,^^^ = 7.1 x 10~^^ for reference. 

4.2.4 BR( Bs ^+^") 

For completeness, we include the bound on the Bs l-i'^ branching ratio [33] 
BR(i?s — > ij,~^iJ:~) < 2.9 X 10^^. It is known that it does not impose any significant 
constraints on the parameter space of mSUGRA. However, for non-universal soft 
terms which we are dealing with, the constraint may be significant [34], especially 
for large tan/3 and low Higgs masses. In practice, we find that the BR(i3s — > /i+^u") 
constraint is satisfied automatically in regions of parameter space allowed by other 
considerations. 

4.3 Example 

Before going into a detailed discussion of our results let us present an example of the 
parameter space allowed by all the constraints. Fig. |51 displays the surviving region 
in the plane (a, m^^2) tan/3 = 5 and /i > 0. The area with a < 5 or so is excluded 
by the presence of tachyons and absence of electroweak symmetry breaking. On the 
other hand, a large a region corresponding to the modulus dominated SUSY breaking 
is excluded by excessive dark matter abundance. The accelerator constraints yield a 
lower bound on the gravitino mass, ~ 30TeV. Very large values of m3/2, except 
perhaps for a very thin strip, are excluded by a combination of the dark matter 
and electroweak symmetry breaking constraints. In Table we provide the SUSY 
spectrum for 3 representative points A,B,C in the surviving parameter space. These 
points are chosen such that the resulting dark matter abundance is consistent with 
the upper and lower WMAP bounds ("favoured" neutralino abundance). 



4.4 Numerical results 

Our numerical results are summarized in Figs. 1 ^111 These plots display contours 
corresponding to various constraints in the (?7i3/2, a) plane for tan (3 = 5, 35, a positive 
and a negative //-parameter and mt = 174, 178, 182 GeV. In addition. Figs. [71 and 
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5466 


4946 


2792 




5902 


5303 


2988 




5761 
5951 


5237 
5453 


2955 
3076 




4662 


3644 


1974 




4669 


3784 


2061 




4332 


3470 


1881 




4507 


3701 


2017 




4506 


3700 


2016 




0.099 


0.105 


0.094 



Table 1: Sample spectra. All masses are in GeV, except for 7713/2 which is given in TeV. 
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tan (3=5 ^>0 mt = 178GeV 




LEP II exclusion 

5 oc 10 

lire 2: Example of the parameter region consistent with the theoretical and accelerator constraints. 



El show the evolution of the SUSY spectrum and the neutrahno relic density with a 
when m3/2 has been fixed. 

There are several features of our analysis that are insensitive to tan/3 and sgn//. 
In all of the considered cases, the parameter space with a ^ 4 is excluded by requiring 
absence of tachyons. As explained in Section 3, a ~ corresponds to pure anomaly 
mediation which predicts tachyonic sleptons. For 4 > a > 2, the sleptons have 
positive masses squared but the squarks turn tachyonic. This feature is specific to our 
scenario and appears due to the mixed modulus-anomaly contribution proportional 
to a in Eq. (ITHl) . 

Another robust feature is the presence of a "no electroweak symmetry breaking" 
(NO EWSB) region adjacent to the tachyonic area. As elaborated in subsection 4.1.1, 
it is related to the suppression of the gluino mass at low a. Electroweak symmetry 
breaking occurs when the Higgs mass parameter m^^^ is negative and sufficiently large 
in magnitude (Eq. (|25|) ). The RG evolution of vr?^^ is controlled to a large extent 
by M3 and for small gluino masses electroweak symmetry breaking is not possible. 
Appearance of the NO EWSB exclusion region at higher 7713/2 when a is fixed is 
associated with two loop effects. 

In what follows, we study effects specific to certain regions of the parameter space. 

4.4.1 Low tan/3 regime 

For small tan/?, the most important accelerator bound is that on the lightest Higgs 
boson mass, see Fig. El (the green (light grey) dashed line). It sets a lower bound on 
the gravitino mass m-3/2 ~ 30TeV which translates into a lower bound on the squark 
and slepton masses of order 2 TeV. 
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Clearly, such a heavy spectrum cannot explain the muon g — 2 anomaly. For 
reference, we display in Fig. Elthe contour (black dashed) corresponding to Ja^^^'^ = 
7.1 X 10^^'^. Above this line the SUSY contribution is too small to be relevant to the 
muon anomaly, yet it is allowed at a 3a level. 

The region between the two solid black contours satisfies the upper and lower 
WMAP bounds. The area above the contours corresponds to excessive neutralino 
abundance and is ruled out, whereas that below the contours is allowed given ad- 
ditional non-SUSY components of dark matter. Fig. [3 explains the shape of the 
allowed region by tracking the composition of the LSP and the SUSY spectrum as a 
function of a at fixed m3/2- For a ~ 5 we are close to the NO EWSB region so that the 
fi term is small and the neutralino LSP is mainly a higgsino. Since m^o ~ "1-^+ ~ 

the coannihilation with the chargino xt work and, furthermore, the higgsino 

coupling to the Z allows for the efficient s-channel annihilation XiXi Z ^ ff. 
This produces acceptable LSP relic abundance. As we increase a, the higgsino gets 
heavier and the LSP becomes more and more bino~like. The annihilation cross sec- 
tion decreases and the relic abundance becomes excessive already at a ~ 5.5. As 
we go to even higher a ~ 7, the neutralino mass approaches the value Ma/2, where 
the annihilation proceeds efficiently through the pseudo-scalar Higgs exchange. As 
a result, the relic density falls and the WMAP bounds are respected. At a > 8 no 
efficient annihilation channel is available and dark matter is overproduced. 

It is interesting to remark here that, in contrast with mSUGRA, the A-pole 
annihilation opens up for a higgsino-like neutralino, and not a 6mo-like one. Similar 
merging between the higgsino and the higgs annihilation branches has recently been 
observed in effective supergravity models with non-universal gaugino masses [36]. 
This effect allows for a large zone of the parameter space respecting the WMAP 
bounds at tan /3 = 5 through efficient annihilation-coannihilation processes, which is 
not the case in mSUGRA. 

4.4.2 Large tan/? regime 

At large tan/3, the Higgs mass constraint becomes less stringent. The main reason 
is that at tree-level ruh ~ Mz\cos2f3\, which increases with tan /?. It requires (Fig. 
ISJ ^-3/2 > 20 — 30 TeV as compared to 7713/2 > 30 — 50 TeV for low tan /3. 

On the other hand, the b ^ sj amplitude grows with tan/3 and becomes more 
important. In particular, it provides the most severe accelerator bound for fi < 
and excludes a large portion of the parameter space. 

We observe that at large tan (3 even a heavy SUSY spectrum contributes signif- 
icantly to the muon g — 2 and ameliorates the anomaly for ^ > (Fig. |H1 region 
below the black dashed line). 

BR(i?s — > iJ-'^fJ,") does not provide any considerable constraint. Indeed, this 
observable is relevant for highly non-universal cases [34] which in our model are 
excluded by the presence of tachyons (low a regime) . 

The evolution of the relic density with a (Fig. ^ differs from the low tan (3 case 
because the pseudo-scalar Higgs exchange followed by a Higgs decay into bb pairs now 
dominates the relic abundance calculation. This process is efficient at large tan f3 for 
two reasons. First, the pseudo-scalar Higgs mass m\ ~ i^'h^ ~ much smaller 

due to the negative bottom quark Yukawa RG contribution to rn^u^. For example, 
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niA ~ 2;U ~ 2m^o already at a ~ 5 (Fig. |H1 left). Second, the Abb coupling is 
proportional to tan j3 and the corresponding cross section cx^^^^^^^ grows as tan^ (5. 
As a result, the relic abundance is well below the WMAP range for 4 < a < 7. As the 
bino component of the neutralino increases, the relic density grows to its maximum 
value around a ~ 9. Then it drops again for a ~ 10 where ttia ~ 2Mi ~ 2m^o 
corresponding to an opening of a bino-like A-pole. 



4.4.3 Influence of the sign of ^ 

It is well known that the b ^ constraint is more important for ^ < [37]. The 
reason is that in this case the SUSY contributions interfere constructively with those 
of the SM increasing the branching ratio, especially at large tan This effect is 
clearly seen from FigsEl and |H1 (right). 

The SUSY contribution to the muon g — 2 usually has the same sign as fi. Thus, 
a positive /i is preferred by the muon anomaly. In any case, the g — 2 discrepancy is 
not conclusive and we treat it as a 3a constraint. 

Other observables are less sensitive to the sign of fi. 



4.4.4 Uncertainties due to the top mass 

We find that some of the results are very sensitive to the precise value of the top 
mass. To take this into account, we provide the exclusion plots for 3 values of the 
top mass: the central value rrit = 178 GeV and the 2a limits nit = 174, 182 GeV'' 

(Figs. ITnlTTT). 

The top mass affects most of all the Higgs mass bound and the relic density. The 

former is sensitive to nit through the one loop correction 6mf^ oc log ( ^ j • For 

a heavier top, a larger portion of the parameter space is allowed by the Higgs mass 
constraint (Figs. ^land^Q right). 

The neutralino relic density is affected by nit mainly in the A-pole region. There 
the neutralino is typically higgsino-like, m^o ~ ^. The value of /z^ ~ ~^'h2 depends 
strongly on mt via the top Yukawa contributions to ni'jj^. For the same reason, the 
pseudo-scalar Higgs mass is sensitive to mt , ni?^ ~ ni'jj_^ — ni^u^ ■ The net result is that 
for larger nit, a broader A-pole region is available (Fig. Ilflj) . This effect disappears at 
large tan /3 in which case the bottom quark Yukawa decreases ni'jj_^ and, consequently, 

We note that at large tan (3 and mt an "mSUGRA-like" A-pole regime becomes 
available. This is seen in Fig. ^2 (right) at a ~ 10 and ~ 25 TeV. The pole 

corresponds to annihilation of bino-like neutralinos through the pseudo-scalar Higgs. 



4.4.5 Summary 

The above analysis shows that there are considerable regions of parameter space 
where the model is consistent with all the constraints. The most restrictive acceler- 



■^The preliminary CDF/DO update [38] yields mt = 172.7 ± 2.9 GeV. The central value mt = 178 GeV 
used in our analysis is within the 2a interval. For smaller mt the Higgs mass bound is more constraining, 
see Figs. UHl HH 
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ator bounds are due to the Higgs mass constraint and BR(6 sj). In some parts 
of parameter space, the muon g — 2 anomaly is amehorated. A positive /x-parameter 
is preferred, whereas both low and high values of tan /3 are allowed. The results are 
sensitive to the top mass such that its higher values lead to larger allowed regions of 
parameter space. 

In most of the cases considered, the resulting SUSY spectrum is rather heavy. 
This can be understood as follows. The Higgs mass constraint yields a lower bound 
on the stop masses of order 1 TeV. Since all the SUSY masses are controlled by 1713/2, 
this bound implies a large 171^/2 and thus a heavy spectrum. This is different from 
the mSUGRA case where the scalar masses, the gaugino masses and the A-term can 
be varied independently. In spite of the heavy spectrum, the degree of fine-tuning 
to get the right EW breaking scale is similar to that of mSUGRA (< 1%) as it is 
mainly sensitive to the 3rd generation scalar masses. 

We note that if we do not insist on the neutralino being the dominant component 
of dark matter, for large tan/? the spectrum is allowed to be lighter, 300 GeV - 1 
TeV. 

5 Conclusions 

In this work, we have studied SUSY phenomenology of the KKLT-type flux com- 
pactification scenario with the MSSM on D7 branes. This setup leads to a specific 

pattern of the soft masses, with modulus and anomaly mediated contributions being 
comparable, and avoids the cosmological gravitino/moduli problems. 

The parameter space includes 3 continuous variables m3/2, a, tan /3 and a discrete 
parameter sgn/x. The resulting SUSY spectrum is non-universal which distinguishes 
the model from mSUGRA and leads to distinct phenomenology. In particular, the 
neutralino LSP is often higgsino-like such that low tan/3 is allowed by dark matter 
considerations, in addition to the usual large tan /3 regime. We find that all exper- 
imental constraints can be satisfied simultaneously in large portions of parameter 
space. Curiously, a ~ 5 required by the shape of the original KKLT lifting potential 
is consistent with the constraints. 

We find that the SUSY spectrum is required to be quite heavy, typically between 
1 and 5 TeV. Although this has certain merits in relation to the CP and flavour 
problems, it may be challenging to discover the superpartncrs at the LHC. Yet, at 
least part of the parameter space with the squarks and gluinos below roughly 3 TcV 
will be explored. We note also that in some cases the charginos may be long lived due 
to their near degeneracy with the LSP (Table 1), which represents a typical anomaly 
mediation signature [39]. 

Finally, it is encouraging that a theoretical model conceived to address the mod- 
uli stabilization problem turned out to have remarkably healthy phcnomcnological 
properties. 
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ERRATUM 



Confusion in various soft term conventions has lead to an unfortunate error in 
our numerical analysis. Derivation of the soft terms in eq. ()12() assumes that the 
Yukawa couplings between gauginos and matter fields contain the i factor: C = 
—iy/2didjKx"'Q^jT""ipi + ^-c. The numerical codes use a different convention (without 
the i factor) which, effectively, amounts to changing the relative sign between the 
gaugino masses and the ^-terms. 

Correcting this error leads to the following modifications: 

1. The A-terms are typically large at the TeV scale. 

2. Some of the parameter space is excluded by the presence of the stop and stau 
LSP. 

3. The neutralino LSP is usually a bino. 

4. The SUSY spectrum can be lighter (below 1 TeV) in certain cases. 

5. The shape of the allowed parameter space somewhat changes (FiglSJ. 

In spite of these changes many qualitative results remain the same. In particular, 
there is a considerable part of the parameter space, e.g. around a ~ 5, which is 
allowed by all the constraints. Typically, the spectrum is heavy (above 1 TeV), see 
Tab 121 although some exceptions can be found. This is enforced by the Higgs mass 
bound and BR(6 — > 57). Acceptable abundance of dark matter can be obtained either 
due to stop-coannihilation or the A-funnel. These results agree with [40], [41]. 

We thank the authors of [40], K. Choi, V. Lowen and H.P. Nilles for important 
communications. 
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Figure 3: Example of the parameter region consistent with the theoretical and accelerator constraints. 
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Table 2: Sample spectra. All masses are in GeV, except for 
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Figure 4; Constraints on the parameter space (see Fig. 3 of the main paper). 
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Figure 5: Constraints on the parameter space (see Fig. 3 of the main paper). 



A: MSSM RG parameters 



In this appendix, we list the MSSM renormahzation group parameters which appear 
in the soft terms formulae ((T^ . The U{1), SU{2) and SU{2>) gauge couplings are 
denoted by Qa- Here U{\) is GUT normalized and related as gi = y^S/S^fy to the 
hypercharge coupling. 

The beta function coefficients ba are defined as 

dga 1 3 



ha9t ■ (A.l) 



91og/u 167r2 
In the MSSM, 

63 = -3 , 62 = 1 , hi= 33/5 . (A.2) 

The anomalous dimension 7j describes the RG dependence of the wave function 
renormahzation Zj, 

d\ogZi _ 1 

51^ " 8^^* ■ ^ "^^ 

In supersymmetric theories, the following general formula holds: 

^, = 2Y^glC^^{Q,)-Y,\y^? ■ (A.4) 

In the second term, the sum runs over all Yukawa couplings yi involving Qi with 
appropriate color factors included. The quadratic Casimir C2{Qi) takes the following 
values: Cf = 4/3 for the SU{3) fundamental or anti-fundamental representation, 
C| = 3/4 for the SU{2) fundamentals, = qf, where qi is the U{1) charge of Qi. 
The anomalous dimensions of the MSSM fields read: 

7Q, = 8/3gl+3/2gl + l/mgl-{y^ + yl)53,, 

-fu, = 8/3gl+8/15gl-2y^63p, 

7D, = 8/3gl + 2/l5gl-2yl63p, 

IL, = 3/2ff| + 3/IO5? - y?<53p , 

7Ep = 6/5c/i - 2y^(53p , 

= 3l2gl^3l\^g\-3yl-yl, 

^H, = 3/2ff| + 3/IO5? - 3y2 . (A.5) 

We have neglected all the Yukawa couplings except for the diagonal ones involving 
the third generation. 

The soft term formulae also involve 7^ = Stt^ ak^^ - From eq. (|A.4|1 . 

7. = 2Y,9Xc^{Q^) - E \y^W^ ■ (A-6) 

Here by- describes the running of the Yukawa couplings, = i^Vibyi- In the 
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MSSM, 

= -8gl + 3/2gl + ll/SO^f - {yX, + ylhy,)5^p , 

ID, = -8gl + 22/25gt-2ylbyj3p, 
jL, = 3/2gl + 99/50gj-yXjsp, 

jE, = l98/25gt-2yXj3p, 

7i/i = 3/2gl + 99/50gt-3yiby^-yX^ , 

jH, = 3/2gl + 99/50gt-3yXt, (A.7) 

where 

by, = 6y^ + yl-16/3gl-3gl-13/15gl, 

by, = yl + Qyl + yl-\&/3gl-3gl-7/lbgl, 

K = 3yl + Ayl-3gl-9/bgl. (A.8) 

Finally, the soft scalar masses contain a mixed anomaly-modulus contribution pro- 
portional to dr^i which appears due to the T-dependence of the gauge couplings. In 
our model, 

{T + T)dTli = -2Y,glci{Qi) + 3Y,\yiW (A.9) 

a i 

such that 



(r + r)9T7Qp = 


-8/3gl 


- 3/2^2 _ l/30gj + 3(y2 + ^2)^3^ 


iT + T)dTlu, = 


-8/3gl 


- 8/15gf + 6y^6sp , 


{T + T)dTjD, = 


-8/3(?i 


- 2/155? + (iylhp , 


{T + T)dTlL, = 


-3/2gj 


- 3/10(7? + 3y%p , 


{T + T)dTlE, = 




+ Qyl^sp , 


{T + T)dTlH^ = 


-3/2ffi 


- 3/lOgl + 9yl + 3yl , 


{T + T)dTjH2 = 


-3/25i 


-3/105? + 9yf . 
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tanP=5 ^>0 mt = 178GeV tanP=5 ^<0 mt = 178GeV 




5 ry 10 5 ,v 10 



Figure 6: Constraints on the parameter space (a, m3/2) of the KKLT scenario at tan/3 — 5, ^ > (left) 
and /i < (right). The region below the light grey (green) dashed line is excluded by the bound on the 
Higgs mass; below the the dotted line - by the bound on the chargino mass; below the solid line - by 
b ^ sj (not shown on the left plot). The area between the black contours satisfies the WMAP constraint 
0.094 < n^oh^ < 0.129, whereas the region above it is excluded due to excessive LSP relic abundance. The 
black dashed line corresponds to a 2(t hmit 5a^^^^ > 7 x 10"^° (left) and a 3a hmit (Ja^^^"^ > -6 x 10"^° 
(right). The lower black dashed line gives a 2(7 upper bound (5a^^^^ < 47 x 10^^°. BR{B ji^ ^i^) 
imposes only a weak constraint and is not shown. 
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Figure 8: As in Fig. Elbut for tan/3 = 35. 
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Figure 9: SUSY spectrum and the relic density as functions of a for tan/3 = 35, > and 1113/2 = 25 
TeV. 
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Figure 10: As in Fig. Elbut for tan/3 = 5, > 0, to* = 174 GcV (left) and 182 GeV (right). 
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Figure 11: As in Fig. Elbut for tan/3 = 35, > 0, m* = 174 GeV (left) and 182 GeV (right). 
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